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i Observations of neutral hydrogen in spiral galaxies reveal a sharp cutoff in the 

radial density profile at some distance from the center. Using 22 galaxies with known 



HI distributions as an example, we discuss the question of whether this effect can be 
associated exclusively with external ionizing radiation, as is commonly assumed. We 



\ show that before the surface density reaches o"hi < 0.5 A^0/pc^(the same for galaxies 

' of different types), it is hard to expect the gas to be fully ionized by background 

Qh, radiation. For two of 13 galaxies with a sharp drop in the HI profile, the "steepening" 
O . 

^ , can actually be caused by ionization. At the same time, for the remaining galaxies, 

' the observed cutoff in the radial HI profile is closer to the center than if it was 
a consequence of ionization by background radiation and, therefore, it should be 

> '. 

\^ . caused by other factors. 

(N 

\^ ' Article is published in Astronomy Letters, 2012, Volume 38, Number 4, Pages 

in ■ 



222-230. 



I. INTRODUCTION 

Neutral hydrogen in galaxies extends to large distances from the center exceeding 
considerably the optical radius i?25- At the same time, the surface density of neutral hy- 
drogen in spiral galaxies drops sharply at some distance from the center. The ionization of 
a neutral gas by ultraviolet background radiation was considered as an explanation of this 
phenomenon (see Bochkarev and Sunyaev 1977; Maloney 1993; Corbelli and Salpeter 1993). 
However, the number of galaxies considered in the cited papers is small; in addition, fairly 
rough estimates of the volume densities for the gaseous and stellar disks were used in them. 
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The elimination of these shortcomings makes this work topicaL 

At the disk periphery, where there is virtually no star formation, the contribution of young 
stars to ionizing radiation is nearly zero, and the contribution of background radiation is 
decisive. According to Haardt and Madau (1996), quasars at high redshifts are mainly the 
sources of ionizing background radiation. Direct observations of background radiation are 
difficult to perform, because it is absorbed by a thick layer of atomic hydrogen in the Galactic 
disk. 

Starting from some distance from the center, the bulk of the hydrogen must be ionized. 
In this paper, we make an attempt to reconstruct the "total" (neutral+ionized) hydrogen 
profiles to estimate the total surface density in the outer regions of spiral galaxies. We took 
the radial HI profiles from the literature and considered them up to distances from the center 
to which the available information allows us to estimatethe volume densities of HI, H2, and 
stars as well as the thicknesses of the gaseous and stellar disks by the method of Narayan 
and Jog (2002). A brief description of the method, references to the radial HI profiles used, 
and the calculated volume densities and thicknesses of the disks are given in Abramova and 
Zasov (2008, 2011). When solving this problem, we took into account the dark halo, which 
plays a prominent role at the periphery of some galaxies. 

In this paper, we reconstructed the total hydrogen profiles for 22 galaxies with radial 
HI profiles known from observations: ten low-surfacebrightness (LSB) galaxies, one normal- 
brightness galaxy with an extended LSB disk (NGC 289), three SO galaxies, and seven 
normal galaxies (including our Galaxy). For comparison, we performed calculations for the 
spiral galaxy NGC 3198, which was considered previously by Maloney (1993) and Maloney 
and Bland-Hawthorn (1999). For this galaxy, we used the radial column density profiles 
and the hydrogen disk half-thickness from Maloney (1993). The list of the galaxies and the 
adopted distances to them are given in the Table [H 

Here, we consider the basic equations solved in the problem, present the results of our 
calculations, and calculate the radial profiles of atomic and ionized gas for a given, monoton- 
ically decreasing total gas profile. In final section "Discussion and Conclusions" we analyze 
the obtained results. 



Table I: Information about the sample galaxies 



Galaxy 


Type 


Mpc 


(1) 


(2) 


(3) 


Malin 1 


LSB 


330 (Moore and Parker 2007) 


Malin 2 (F568-6) 


LSB 


186 (Pickering et al. 1997) 


F561-1 


LSB 


62.7 (de Block et al. 1996) 


F574-1 


LSB 


96 (Swaters et al. 2003) 


F568-1 


LSB 


85 (Swaters et al. 2003) 


F568-V1 


LSB 


80 (Swaters et al. 2003) 


F568-3 


LSB 


77 (de Block et al. 2001) 


UGC 128 


LSB 


60 (van der Hulst et al. 1993) 


UGC 1230 


LSB 


51 (de Block and Bosnia 2002) 


UGC 6614 


LSB 


85 (de Block et al. 2001) 


NGC 289 


SBbc+LSB 


21 (Walsh et al. 1997) 


UGC 2487 (NGC 1167) 


SO 


67.4 (Noordermeer 2006) 


UGC 11670 (NGC 7013) 


SO 


12.7(Noordermeer 2006) 


UGC 11914 (NGC 7217) 


SO 


14.9 (Noordermeer 2006) 


M33 (NGC 598) 


Sc 


0.7 (Boissier et al. 2004) 


M51 (NGC 5194) 


Sbc 


8.4 (Boissier et al. 2004) 


M81 (NGC 3031) 


Sab 


3.6 (Boissier et al. 2004) 


MlOO (NGC 4321) 


SABb 


17.0 (Boissier et al. 2004) 


MIDI (NGC 5457) 


SABc 


7.5 (Boissier ct al. 2004) 


M106 (NGC 4258) 


SABb 


8.0 (Boissier et al. 2004) 


Galaxy 


SBbc 




NGC 3198 


SB(rs)c 


9.4 (Maloney 1993) 
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II. BASIC EQUATIONS 

The main adopted simplifications are as follows. The HI disk is assumed to be an 
equilibrium and axisymmetric one, and the thickness of the stellar disk is determined by the 
condition of gravitational stability that imposes constraints on the stellar velocity dispersion. 
In addition, we disregard the local gas density inhomogeneities, whose presence changes the 
gas ionization fraction. 

The ionization balance equation for the gas layer at distance z from the disk midplane 
{z = 0) is 

where z is the coordinate along the axis perpendicular to the disk plane; hv^ = 13.6 eV; 
Ue and Hp [cm~^] are the ion and proton number densities, respectively; nni [cm~^] is the 
number density of hydrogen atoms; — 7 ■ 10~^^ (^^^)^ [cm^] is the photoionization cross 
section for the hydrogen atom; Ju{z) [ cm^^sHzsr ] ionizing radiation at given 

z; aB is the radiative recombination coefficient, which generally depends on the amount of 
neutral hydrogen and the local electron temperature. Since the contribution of elements 
heavier than hydrogen is small in our problem, below we assume that 7?/g ~ Tip. 
According to Ferland (1980), the total recombination coefficient 

oo 
n=2 

can be approximated to within 10% by the expression 



2.9 • 10-1° T^o-^^ Te<2.6-10^K; 
1.31 ■ 10^^ T-i-i^ Te > 2.6 • 10^ K. 



In our case, Te = 10^ K, so that aB = 2.4 • lO^^^ cm^s^^ 

The initial flux of ionizing radiation J° (uncorrected for absorption) can be written as 
(see Bland-Hawthorn et al. 1997) 



erg 



(2) 



-cm^ sHz sr. 

where J°2i = 0.08 t.in2^gHzsr metagalactic flux near the Lyman limit i/ — uq. If 

the cosmic radiation is determined by quasars at high redshifts, then, according to Bland- 
Hawthorn et al. (1997), p^l, i.e., J° = 0.08 • lO-^i {^f) 



cm^ s Hz sr ' 
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When penetrating into the disk, part of the initial flux is absorbed and the attenuated 
flux at distance z from the disk midplane will be 

We will assume that the number densities rie = rip and whi decrease along the z coordinate 
exponentially: 

rip,ni = nl^je-^, (4) 

where the vertical scale height of the disk Zg is related to the half-thicknesses of the disks of 
the sample galaxies hi/2 [pc] by the relation 

Zs = 3.086 • 10^^ (5) 

The values of hi/2{z) for the galaxies considered here were calculated previously (Abramova 
and Zasov 2008, 2011). 

For the convenience of our calculations, we will change the system of units and pass from 
the particle number densities ni to the column densities Ni [cm~^], which are related to the 
observed quantities, the surface density cij, by the relation 

N, = 1.26 X lO^Vi, where (6) 

(Tj are measured in A^o/pc^. Let us express Ni in terms of the particle number densities rii 
(the subscript i can refer to both protons p and neutral HI atoms): 



OO POO 



Ni= ni{z)dz = 2 ni{z) dz = 271^^^ e"~ ciz = 2n°z„ (7) 

J-oo Jo Jo 

where are the number densities of the corresponding particles in the disk midplane at 
z = 0. Thus, if Hp and nni change exponentially along the z coordinate, then Ni is expressed 
linearly in terms of n^. Since the ionization balance equation ([1]) is valid for all z, we will 
consider it for the case of 2; = 0. 

A. Equation ([T]) for the Disk Midplane 

For z = 0, the flux of ionizing radiation ([3]) attenuated by absorption can be written 

as 

Mz) = J,{0) = 2jy^'^<^''. (8) 
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The factor 2 on the right-hand side of ([8]) means that the flux comes from two directions. 
In the ionization balance equation ([1]), we will pass from integration over the frequency to 
integration over the energy E = hu and from the particle number densities rii to the proton 
(ionized gas) and HI surface densities. Given (E]), ([5]), and ([8]), we will write ([T]) for dp: 

i.i.io«9? 

a^ = 2.1-10VHi/ii/2 / dE. (9) 

J13.6 

Let us denote A = 1.1 ■ IO^cthi and Eq = 13.6 eV. The integral on the right-hand side 
of (El) will then be rewritten as 



^ A A 

_ 1 



dE=^, '</z'e-^ dz'. (10) 



According to Gradshtein and Ryzhik (1963), 

^u-l g-M^ ^ ^-v ^^^.^ ^^-j^ 



where 7(z/; /iu) is the incomplete gamma function. 

The final ionization balance equation in the disk midplane written for the sought-for 
quantity, the proton (ionized gas) surface density, will take the form 



.,^ = 7.10^f2l|^,(l;^).^=1.1.10».„,. (11) 



The right-hand side of (jTTj) depends only on two quantities: observed surface density of 
atomic hydrogen cthi and the half-thicknesses of the gaseous disks of the galaxies considered 
hi which were found for all galaxies of the sample previously (Abramova and Zasov 2008, 
2011). 

Equation (fTT]) is solved by numerical methods. 



III. THE RADIAL PROFILES OF TOTAL (NEUTRAL + IONIZED) AND 

IONIZED HYDROGEN 

Figures ([HIS]) present the modeled radial surface density profiles for ionized and total 
(neutral + ionized) gas. The objects of all types belong to the galaxies that have a cutoff in 
the observed radial HI profile: LSB (F561-1, F568-vl, F568-6, F574-1, UGC 128, UGC 1230, 
UGC 6614, Malin 1), SO (UGC 2487, UGC 11914), and normal galaxies (M33, MlOO, and 
NGC 3198). For M81 there may be a cutoff in the radial HI profile at the boundary of the 
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Figure 1: Modeling of the radial surface density distribution for ionized and total (ionized+HI) gas based 
on the observed radial HI surface density distribution for LSB galaxies: the solid line, black asterisks, and 
gray circles indicate, respectively, the total gas density, the ionized gas distribution, and the radial HI 

distribution. 

region considered here, while for the remaining galaxies a gradual decrease in the amount 
of atomic hydrogen with increasing distance from the center is predominantly observed. 

We see from Figs, (^j^ that the amounts of ionized and neutral gas for the galaxies of 
all types become equal in the range of surface densities cthi = 0.3 ^ 0.5 A^o/pc^ (utot = 
0.6 ^ 1.0 A^0/pc^), i.e., when the column density reaches A^hi = 3.8 • 10^^ ^ 6.3 • 10^^ cm~^. 
Both among the models proposed by Maloney (1993) for NGC 3198 and among the models 
computed by Corbelli and Salpeter (1993) for M33, there are those that are consistent with 
this result. 

For the normal galaxies (Fig. [2]) in the R range considered here, the relative ionized gas 
fraction is small. According to Narayan and Jog (2002), the amount of HI in the range 
R = 3 — 16 kpc in our Galaxy (Fig. [2]) is almost constant and the ionization grows weakly. 
The behavior of the giant normal+LSB galaxy NGC 289 is identical to that of the normal 
galaxies considered. The SO galaxies UGC 2487, UGC 11670, and UGC 11914 (Fig. ED, as 
in the case of densities (see Abramova and Zasov 2011), are closer in their properties to the 
LSB galaxies. 

In the case of LSB and SO galaxies, for which the gas ionization is more significant (see 
Figs. ([1]) and ([3])), and in the normal galaxy MlOO (see Fig. |2]), the ionization fraction near 
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Figure 2: Same as Fig. ((T]) for normal galaxies and the normal+LSB galaxy NGC 289. 




Figure 3: Same as Figs. Q and ^ for SO galaxies (UGC 2487, UGC 11670 and UGC 11914), two giant 
LSB galaxies, Malin 1 and F568-6 (Malin 2), and the spiral galaxy NGC 3198 considered previously by 
Maloney (1993) and Maloney and Bland- Hawthorn (1999). 

the cutoff in the observed radial HI profile is low. This may suggest that the sharp drop in 
the observed HI density is caused not by ionization but by a sharp drop in the initial gas 
density. 



9 



IV. CALCULATION OF THE RADIAL HI AND IONIZED GAS PROFILES 
FROM THE KNOWN TOTAL GAS DENSITY PROFILE 

A. The Model Galaxy 

Suppose that the initial density profile is known. Let us write the ionized gas surface 
density as dp = o"tot — <^m- We will assume the total gas density (cTp + am) to be known. 
We solve (fTT]) numerically for cthi and use the value found to calculate the radial ionized gas 
surface density profile. 

Consider a model galaxy whose parameters are chosen arbitrarily based on typical values 
for the normal galaxies considered here. Let the total gas surface density in the outer region 
of the model galaxy drop exponentially 

(Ttot, [-Mo/pc'] = lOe^, R > 5 kpc, (12) 

and the half-thickness of the gaseous disk increase linearly with radius 

hi/2{R), [pc] = 90R - 250. (13) 

Figure (jl]) presents the results of our calculations. The dashed lines indicate the column 
densities at which (according to the estimates from Abramova and Zasov (2011)) hydrogen 
turns out to be ionized over the entire thickness of the gas layer (the socalled critical density). 
We see that in the range of critical densities the mass fraction of ionized hydrogen is 30-50%. 

The cutoff in the theoretical HI profile begins where the densities of atomic and ionized 
hydrogen in the model galaxy become equal. This occurs near cthi ~ 0.5 A^q/pc^, which 
is in good agreement with the results obtained for the 22 galaxies considered here. When 
the amounts of atomic and ionized hydrogen become equal, a transition zone with an extent 
of several kpc in which the solution is mathematically unstable begins (hatched in Fig. (jll). 
In the transition zone, the gas ionization fraction increases sharply with distance from the 
center (see Fig. (jlj)). 

B. Examples of Modeling Several Galaxies 

As an illustration, let us solve the inverse problem considered above for several galaxies 
with farreaching am profiles. Let us take the two LSB galaxies UGC 128 and UGC 1230, 
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R, kpc 

Figure 4: Modeling of the radial atomic and ionized hydrogen profiles for an exponential initial gas 
surface density distribution. The dashed lines indicate the critical densities; the hatching marks the 
regions of mathematical instability of the solution (see the text). 

the giant lenticular galaxy UGC 2487, and the spiral galaxy NGC 3198. We will specify the 
radial total gas profiles in the form of an exponential with a radial scale length corresponding 
to the observed profile in the inner region of the galaxy and will solve Eq. (ITT]) for HI. 

For each of the four galaxies, Fig. ([5]) shows the model exponential initial gas profiles as 
well as the observed and theoretical atomic hydrogen profiles. 

We see from Fig. ([5]) that for two galaxies (the giant SO UGC 2487 and the spiral 
NGC 3198), when the densities reach their critical values Nc, the HI densities change 
abruptly, just as in the model galaxy considered in the previous section. This occurs near 
(Thi ~ 0.4 Ai^/pc^, which is also in good agreement with the results obtained for the 22 
galaxies considered here and with the results obtained for the model galaxy. 

For three galaxies (UGC 128, UGC 1230, and UGC 2487), the HI cutoff profiles do not 
correspond to that expected under ionization: for the LSB galaxies UGC 128 and UGC 1230, 
which are characterized by higher HI surface densities, the model curves of the radial HI 
distribution lie on the graphs above the observed ones and, in contrast to the observed 
curves, have no distinct cutoff, while the theoretical HI profile for the giant SO UGC 2487 
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Figure 5: Reconstruction of the radial HI profile from the specified exponential total gas profile for the 
four galaxies. The dashed lines indicate the critical densities; the regions of mathematical instability of the 

solution are shaded (see the text). 



lies below the observed one and, in contrast to it, has a jump. The result of modeling 
NGC 3198 agrees well with the observations. 

Thus, in the galaxies where the HI surface densities at the edge of the optical disk reach 
o"Hi < 0.5 A^o/pc^, the sharp cutoff in the observed HI profile can be caused by the ionization 
of gas by background radiation. At the same time, for galaxies with higher surface densities 
at the disk edge, the steepening of the observed HI profile begins closer to the center than 
could be expected in the case of ionization by background radiation and, therefore, must be 
caused by other factors. Note that Portas et al. (2009), who analyzed the distribution of 
atomic gas in nine galaxies, made a similar suggestion: the radial HI profiles are probably 
determined by internal causes rather than external ones, to which the background radiation 
belongs. 
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V. DISCUSSION AND CONCLUSIONS 

When solving the problem, we made a number of simplifications. We assumed that 
there were no inhomogeneities in the gaseous disk (allowance for them reduces the ioniza- 
tion fraction). In addition, we disregarded the fact that for the regions deep in the disk, 
the ionizing radiation loses its isotropy. These simplifications could lead the effect that is 
observed in the galaxies Malin 1 and NGC 3198 (Fig. |3]), for which the total gas profile at 
the disk edge turns out to be bent upward. This bend is unlikely to be real; this is most 
likely the result of an overestimation of the ionized gas mass at the edges of the disks of 
these galaxies. 

The overestimation of the ionized gas density can be not only the result of the above 
simplifications but also, in the case of Malin 1, a consequence of the peculiarities of the 
numerical method used to determine the gaseous disk thickness (see Abramova and Zasov 
2011). Indeed, in the regions of very low densities < 10"^^ g/cm^ characteristic of the edge of 
the Malin 1 disk, the accuracy of the numerical method used is insufficient. This could lead 
to an overestimation of /ii/2 at the disk edge and, as a consequence, to an overestimation 
of the ionization and a bend of the theoretical total gas density profile. The possibility 
that the half-thickness of the gaseous disk of NGC 3198 taken from Maloney (1993) was 
overestimated must not be ruled out either. This could also be responsible for bend of the 
theoretical total gas density profile in this galaxy. 

The results of solving the direct and inverse problems for galaxies of different types 
(LSB, SO, and normal) lead us to the following conclusion: for galaxies in which the sharp 
steepening of the radial HI profile begins in regions with densities cthi < 0.5 A^o/pc^, it 
can actually be caused by the ionization of gas at the disk edge. In this paper, two of the 
galaxies considered here, the giant SO UGC 2487 and the spiral NGC 3198, belong to such 
objects. At the same time, for the remaining galaxies with a cutoff in the radial HI profile, 
the sharp decrease in the observed amount of atomic hydrogen in the disk begins well before 
the gas surface density reaches cthi ~ 0.5 A^o/pc^ and, consequently, the cutoff in the radial 
HI profile in them is not related to the presence of ionizing radiation but is attributable to 
other factors. 

Even a threefold increase in the fiux of ionizing radiation does not change the situation 
(see Fig. [6]). This figure presents the solution of the direct problem for four galaxies of 
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Figure 6: Results of modeling the radial surface density distribution for ionized and total (ionized+HI) 
gas based on the observed radial HI surface density distribution for a threefold increase in the flux of 
ionizing radiation The notation is the same as that in Figs. ([TJ31)- 



different types. The increase in the radiation ffux had no noticeable effect on the ionization 
of gas in the normal galaxy MlOO, the normal+LSB galaxy NGC 289, and the SO galaxy 
UGC 11914. At the same time, the ionized gas fraction in the LSB galaxy F568-vl increased 
to an extent that its amount at the disk edge exceeded the amount of HI (cf. Fig. [1]). 
However, despite this fact, the result qualitatively remained the same for all four galaxies: 
the steepening of the observed HI profile begins closer to the center than might be expected 
in the case of background ionization. 

The inhomogeneity of the interstellar medium is an important factor that can change the 
estimate of the gas ionization fraction in disk regions far from the center. The concentration 
of gas into clouds or filaments will reduce its ionization fraction, moving the HI profile cutoff 
even farther from the center, and cannot explain the sharp decrease in density that begins 
earlier than that expected for themean density. However, this effect is difficult to estimate 
quantitatively because of the lack of sufficient information about the small-scale distribution 
of gas. 
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Basically, the increase in ionization fraction could be explained by internal sources of 
ionizing radiation (young stars, supernovae), but the available data for the galaxies under 
consideration give no grounds for this. Indeed, most of the galaxies with an observed steep- 
ening of the HI profile belong to lenticular or LSB galaxies, and, consequently, they have 
neither contrasting spiral arms nor intense star formation. Only three of these galaxies 
(M33, MlOO, and NGC 3198) belong to late morphological types; for the latter, no internal 
sources of ionization are required to explain the cutoff in the HI profile. In addition, the 
spiral structure at the periphery of galaxies (in the regions where the cutoffs are observed) 
is rather indistinct. 

Note also that among the 22 galaxies considered, there are no objects located in clusters 
(except for MlOO at the periphery of the Virgo cluster) and only four galaxies (M81, M33, 
NGC 289, and UGC 2487) are members of groups. No cutoff in the HI distribution is 
observed in the range of radii considered in the galaxy with a close companion (M51). 

Therefore, we cannot consider the ionization of a rarefied gas at the edges of galactic 
disks by external radiation as a universal mechanism leading to the observed cutoff in the 
radial HI profiles, and other explanations for this phenomenon related to the formation and 
evolution of disks should be sought. 
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